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^ : ABSTRACT 

^ ' Although analytical studies on the secular motion of the irregular satellites 

>■ ■ have been published recently (Kinoshita & Nakai 1999; Yokoyama et al. 2003), 

these theories have not yet been satisfactorily reconciled with the results of di- 
^ ' rect numerical integrations (Carruba et al. 2002; Nesvorny et al. 2003). These 

Q ■ discrepancies occur because in secular theories the disturbing function is gen- 

^ ■ erally averaged over the Sun's orbital motion, whereas instead one should take 

' into account some periodic terms, most notably the so-called "evection", which 

P^' can be large for distant, slow- moving satellites. This problem is identical to that 

Q ■ initially encountered by Newton and other historical researchers when studying 

-t— > ■ the Moon's motion. Here we demonstrate that the evection and other terms 

a • from lunar theory can be incorporated into the more modern Kozai formalism, 

^ ■ and that our synthetic approach produces much better agreement with results 

^ ■ from symplectic integrations. Using this method, we plot the locations of secular 

■ resonances in the orbital-element space inhabited by the irregular satellites. Our 

model is found to predict correctly those satellites that are resonant or near- 
resonant. 

We also analyze the octupole term in the disturbing function (Yokoyama et 
al. 2003) to determine the strengths of resonant-locking for satellites whose lon- 
gitudes of pericenter are librating. By independently integrating these satellites' 
nominal orbits using a symplectic integrator, we show that the strength of this 
resonance can be successfully obtained from simple analytical arguments. 

We note that the distribution of irregular satellite clusters in the space of 
proper orbital elements appears to be non-random. We find that the large ma- 
jority of irregular-satellite groups cluster close to the secular resonances, with 
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several objects (Pasiphae, Sinope, Siarnaq, formerly S/2000 S3, and Stephano, 
formerly S/1999 U2) having practically stationary pericenters. After proposing 
the name "Main Sequence" to describe this grouping, we point out that none of 
the largest satellites (those with radii i? > 100 km) belong to this class. Finally, 
we argue that this dichotomy implies that the smaller near-resonant satellites 
might have been captured differently than the largest irregulars. 

Subject headings: celestial mechanics — planets and satellites: general — planets 
and satellites: individual(Pasiphae, Siarnaq, Sinope, Stephano) 

1. Introduction 

Irregular satellites are usually considered to be those that orbit the giant planets far 
beyond the major moons, often with quite eccentric and inclined orbits. The number of 
known irregular satellites has grown rapidly (Gladman et al. 1998, 2000, 2001; Sheppard 
& Jewitt 2003; Holman et al. 2004) at the turn of the 21^* century due to advances in 
CCDs and computing, coupled with the availability of numerous large telescopes. The new 
objects tend to group into clusters (Gladman et al. 2001) that arc, at least in the case 
of the Jovian satellites, centered on large objects that have been known for most of the 
last century. While modern observers have improved upon the efforts of their predecessors 
(and sometimes spectacularly so, with the year 2000 re-discovery of Themisto, that was lost 
since 1975), the theoreticians have not been similarly active in bettering analytical models, 
instead generally preferring numerical simulations for determining histories (Carruba et al. 
2002, 2004; Nesvorny et al. 2003; Astakhov et al. 2003). Thus, for the most part, papers 
from two generations ago (Brown 1930; Brown & Brouwer 1937) which addressed Pasiphae's 
motion analytically, have been almost forgotten, even though they treat this problem in 
more detail than any of the modern work. Of all the contemporary researchers, only Saha 
& Tremaine (1993), while concentrating on a numerical approach, have correctly attributed 
Pasiphae's secular^ resonance to mutual cancellation between various terms that appear 
in the expression for the precession of apsides that was obtained in lunar theory. None 
of the more recent papers (Kinoshita & Nakai 1999; Carruba et al. 2002; Yokoyama et 
al. 2003; Nesvorny et al. 2003) reference Brown's work, nor do they mention evection and 
other corrections to secular behavior that arise when the disturbing function is appropriately 
averaged. 



^The term "secular" refers to averaged or non-periodic behavior, and more specifically to perturbations 
which occur on a precessional timescale rather than on one comparable to an orbital period. 
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The difficulty in describing the secular behavior of a distant sateUite is most famously 
illustrated by Newton's inability to explain the precession of the lunar apsides. The pre- 
cession rate of an ecliptic, near-circular satellite orbit under the solar influence, to the first 
order, is given by 



where w is the satellite's longitude of pericenter, Up is the planet's mean motion, Ug is that 
of the satellite and the dot signifies differentiation with respect to time (Danby 1992). For 
Earth's Moon, this formula gives an apsidal precession period of about 18 yrs, identical 
to that of the nodal regression. However, the observed precession of the lunar apsides is 
twice as fast. Newton himself was frustrated by this discrepancy, and considered it to be 
a significant faihng for his theory of gravitation. Clairaut, in 1847, reahzed that this large 
error in precession rate arises from evection terms which remain after averaging over the 
orbital motion has been completed (Baum & Sheehan 1997). 

Lunar evection, known to Hipparchus, and nicely treated by Brouwer & Clemence 
(1961), is, in essence, a periodic perturbation to a satellite's orbit that depends on the 
angle between the satellite's fine of apsides and the Sun's instantaneous position. The rela- 
tive strength of this perturbation is determined by m = rip /rig, the ratio of the mean motions 
of the planet and the satellite (see Eq. 1). For Earth's Moon, m is about 1/13, while it can 
be as large as 1/5 for some retrograde Jovian satellites. In the next section, we show that, 
for an m of only 0.1, evection-induced precession of the apsides surpasses in magnitude that 
arising from purely secular terms. 

When the Space Age started nearly fifty years ago, the focus for the hierarchical three- 
body problem (Sun, planet, satellite) veered away from the motions of the Moon and other 
natural satellites to the orbits of artificial satellites; more recently, planets in binary systems 
(Innanen et al. 1997; Holman & Wiegcrt 1999) were studied. This shift has redirected 
the emphasis from a very accurate description of orbits that are close to being circular 
and equatorial, to a more qualitative treatment of orbits having arbitrary inclinations and 
eccentricities. Generally, numerical integrators have been employed whenever the celestial 
bodies' positions need to be accurately predicted. 

Analytical models of the hierarchical three-body problem usually average the leading 
term (in tts/ttp ~ m^/^) of the disturbing function over the orbital motions of both perturber 
and the perturbee, and obtain the differential equations for the evolutions of inclination i, 
eccentricity e, argument of pericenter uj and longitude of ascending node Q for arbitrary 
eccentricity and inclination (Innanen et al. 1997, corrected by Carruba et al. 2003): 
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^ = — ^e^sm(2a;) sini cosi/Vl — e^, (2) 



de 15 



, = — e Vl — e^sin(2a;)sin^i, (3) 
or 8 



^ = ^{2(1 - e^) + 5 sm''u;[e^ - sm^y/VT^, (4) 
or 4 



and 



^ ^ _^j3 ^ ;^2e2 - 15e2 cos^u]/VT^, (5) 
dr 4 

where time has been non-dimensionalized by introducing t — {1 — ep)~^^^mnpt, with Cp 
being the planet's heliocentric eccentricity. The above equations can either be integrated 
numerically (Carruba et al. 2002; Ncsvorny ct al. 2003) or their solutions can be expressed 
in terms of elliptical integrals (Kinoshita & Nakai 1999). 

A satellite's secular behavior obtained from this low-order theory can be divided into 
cases when the argument of pericenter circulates or librates. According to this simplified 
theory, for all satellites with inclinations below 39.2° (or greater than 140.8°), u circulates 
from to 27r, while for orbits with higher inclinations (39.2° < i < 140.8°) lv either cir- 
culates or librates around 90° or 270°. The last case is what is usually referred to as the 
"Kozai resonance" , and was among the most important discoveries of the "new celestial me- 
chanics" for the hierarchical three-body problem (Kozai 1962). The first natural satellites 
observed to be in the Kozai resonance were Saturn's Kiviuq and Ijiraq (Vashov'yak 2001; 
Cuk et al. 2002a; Nesvorny et al. 2003; Yokoyama et al. 2003); now about a half-dozen are 
known. 

However, as we have gradually come to realize, the Kozai theory, as developed so far, 
is seriously limited. We have already mentioned that, for small eccentricity and inclination, 
the formula for the precession of w (defined as Q + u; for prograde bodies) reduces to Eq. 
1, and gives an error of 50% for the Moon. Furthermore, the Kozai model suggests that 
prograde and retrograde satellites should have symmetric secular behavior: yet, while the 
Moon's apsides precess faster than they "should" , retrograde satellites like Jupiter's Pasiphae 
and Sinope show very slow precessions of their apsides, and it is this that enables a secular 
resonance to take effect (Whipple & Shelus 1993; Saha & Tremaine 1993). Even though 
the strength of the locking between the planet's and the satellite's lines of apsides depends 
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on the octupole term in the disturbing function (Yokoyama et al. 2003; Sec. 4 below), the 
resonance first requires that the leading terms give tzj ~ 0, which is impossible according to 
Eq. 1 and becomes possible only for a specific inclination (~ 36° for e ~ 0) in Eqs. 2 to 5. 
It is historically interesting to note that Yoshihide Kozai himself, in the very last paragraph 
of his landmark paper (Kozai 1962), mentioned that his method cannot be applied to the 
most distant satellites, due to problems with averaging over the solar true longitude. 

Many researchers have worked on the lunar theory since Newton and Clairaut (Tisserand 
1894; Brown 1896; Brouwcr & Clemence 1961), and - for this case of low e and i - they have 
obtained an expression for ikj that is a high-order polynomial in m (Tisserand 1894): 



This expression makes clear that there may be (negative) m's for which w = 0. Saha 
& Tremaine (1993, their Fig. 1) recognized Eq. 6's relevance in determining the secular 
resonances of Pasiphae and Sinope. The second term, which invariably produces prograde 
motion, derives solely from evection, as we will show next. Although evection is the largest 
correction to the purely secular motion, the convergence of the series is slow (because higher- 
power terms in m have larger coefficients as seen in Eq. 6) and usually terms with powers 
of m as high as 10 are listed. Higher-order terms come from different sources (mostly from 
periodic terms that do not completely average out) and their calculation requires a quite 
specialized theory. 

Unfortunately, the classical approaches used to derive Eq. 6 are not practical for 
irregular-satellite applications. All lunar theories assume small eccentricity and inclination, 
and arc always constructed in specific variables which arc significantly removed from the 
usual orbital elements. The only instances when the "old celestial mechanics" was applied 
to the study of an irregular satellite's orbit were a pair of papers by E.W. Brown on Pasiphae 
(1930, Brown & Brouwer 1937). In our opinion. Brown's approach is remarkably difficult, 
yet yields only approximate results. 

The main goal of this paper is to incorporate evection and other terms from Eq. 6 into 
Kozai's set of equations. To do this, in the next chapter we explicitly derive the dependence 
of the evection term on eccentricity and inclination, and then argue how other, higher-order 
terms should behave with inclination. In the subsequent chapters we compare our results 
to those from direct numerical integrations, and discuss the possible significance of secular 
resonances for the origin of the irregular satellites. 





-6- 



2. Development of the Disturbing Function 

As usual, we consider the planet to be the central body and the Sun to be an exterior 
perturber. In that case, the disturbing function describing the solar influence on the satellite 
is (Murray & Dermott 1999): 

' °° I 

1=2 

where fj,' is solar mass times the gravitational constant, r and r' are the distances of the 
satellite and the Sun, respectively, from the planet, while ip is the angle between their 
position vectors. Pi{cosip) are Legendre polynomials of order / with the argument cos?/'. 
Since the ratio r/r' is small, of 0(m^/^), for most applications only the leading (Z = 2, or 
quadrupole) term in the disturbing function is taken. The next term (/ = 3, or octupole), 
which is important only within secular resonances, will be discussed later. The quadrupole 
term is 

^ = l^'{?^) ^(3cosV-l)- (8) 
Prom spherical trigonometry, we can expand the angle ip as: 

cos'^ = cos(A' — il) cos(a; + /) + sin(A' — O) sin(a; + /) cosi, (9) 

where / is the satellite's true anomaly, and its inclination i is measured from planet's orbit 
plane. Substituting Eq. (9) into (8), 

R^H'(^L-'^ 1[3 cos2(A' - n) cos^{uj + /) + 3 sin2(A' - n) siii^ {uj + f) cos^ i 
+6 cos(A' — fl) cos{u + /) sin(A' — fl) sin(a; + /) cosi — 1]. 

For practical reasons, we divide this expression into distinct terms: 

u' 1 

R^-^ -(3Ri + 3R2Cos^i + 6Ri2Cosi- Ro), (10) 
r' 2 

where 

R^ = cos^(A' - fl) [cos^ (J cos^ / - | sin(2a;) sin(2/) + sin^ u; sin^ /], 
R2 = sin2(A' - fl) [cos^ u; sin^ / + | sin(2u;) sin(2/) + sin^ u; cos^ /], 
i?i2 = isin(2A' - 2fl)[cos{2uj) sin(2/) + sin(2a;) cos(2/)]. 
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and 

Ro = r\ 

Here we have separated the dependence on uo and /, in order to prepare the expression for 
averaging over /. 

We average these expressions over the sateUite's orbital motion by integration: < Ri >= 
^ Ri Usdt. This is most easily done by expressing r, / and Ugdi in terms of the eccentric 
anomaly E. Thus 

< >= 0^(1 + le^) 
< r^cosV >= a?{\ + 2e^) 
<r^sm''f>^ a\l-e^)/2, 

and 

< sin(2/) >= 0. 
Putting the averaged values back into the expressions for Ri: 

<Ri a2cos2(A'-0)[i(l-e2) + | e^cos^cu], 
<R2 >=a2sin2(A'-l])[|(l-e2) + | sin^ to], 
< Ri2 >= isin(2A' - 20) | e2sin(2u;), 

and g 

<Ro >= 0^(1 + -e'). 

Since both the Kozai theory (Eqs. 2 to 5) and the lunar evection (Brouwer & Clemence 
1961) depend on double angles, we express R in terms of double angles only. With these 
average values, Eq. (10) can be rewritten as 

R' = R{2 r'V/xV) 

= I - 3e2 - 1(1 - e2) sin^ i -9.)^ f e^j [1 + cos(2A' - 29)] [1 + cos(2a;)] 

+ cos^ i[l - cos(2A' - 20)] [1 - cos(2cu)] + 2 cosi sin(2c;) sin(2A' - 2Vt)} 

= I - 3e2 - 1(1 - e") sin^ i[l - cos(2A' - 29)] 
+ f e2{(l + cos2i)[l + cos(2cj) cos(2A' - 29)] 
+ (1 - cos^ i)[cos(2A' - 20) + cos(2u;)] + 2 cosi sin(2u;) sin(2A' - 29)}. 

The resulting disturbing function can be divided into three pieces: 

R^R' ^^Rk (2a;) + Rj (2A' - 29) + Re {2uj, 2X' - 29) (11) 
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with the parts being 



Rj, = ^{1 _ - ^(1 - e^) sin^ i + — e'[l + cos^ i + sin^ i cos(2cu)]}, (12) 
= ^£(l-e2) + ^e2]sin2zcos(2A'-2f]), (13) 



-Re = — ^e^{(l + cos^ i) cos(2a;) cos(2A' - 2Q) 



and 

15 ii'a^ 

16 r' 

+2cosisin(2co') sin(2A' - 2Q)}. (14) 

The first term, Rk-, is the only one that is retained in the Kozai theory. Since it is independent 
of A', the averaging over the "orbital motion" of the Sun (which is, of course, only the 
reflection of the planet's motion) is trivial. The task reduces to: 

1 r^-npdt_ 1 ^^^^ 



271 Jo r'3 a'3(l-e'2)3/2- 



When wc put this result back into Eq. 12, introduce b' = a'yl — ef^ (semiminor axis of the 
planet's orbit) and group the terms differently, we get: 



i?^ = ^[2 + 3e2-(3 + 12e'-15e2cos2a;)sin2i]. (16) 

So 

This is the exact form of the disturbing function obtained by Innanen et al. (1997), which, 
through Lagrange's equations, yields Eqs. (2-5). 



Now we consider Re, which, if i = (i.e., cosi = 1), reduces to 

15 yu'a^ 
"8 r 



Re = — cos(2A' - 20 - 2lu) 

^3 e^cos(2A'-2G7); (17) 



15 fi'a^ 



8 r 

this is equivalent to the evection term given by Brouwer & Clemence (1961) in their chapter 
on lunar theory. However, since we are interested in irregular satellite dynamics, we cannot 
restrict i to be 0° and must retain the dependence on i. Also, wc cannot use zu at this 
point, since it is not defined similarly for all satellites: zu = fl + uj for prograde bodies, but 
zu — fl — cu foT retrograde ones. Nevertheless, we can calculate the secular trend resulting 
from Re in much the same way that Brouwer & Clemence (1961) did. We recognize that 
this term produces a periodic perturbation in e. Using the Lagrange equation for e (Danby 
1992), we find 

V dv / E o jJiT 
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-2cosicos(2a;)sm(2A' - 2Q)]. (18) 

If we consider that A' ~ npt, and that it changes much faster than the satelhte's orbital 
elements, we can integrate Eq. 18, getting 

1 ' 3 

Sce = — eVl - e'm + cos^ i) sin(2c.-') sin(2A' - 20) 

16 jinpT' 

+2cosicos(2u;)cos(2A'-2Q)]. (19) 

Since a^//i = n^^ and /x'/r'^ ~ n^, the factor [(/x'//x)(ns/np)(a/r')^] in Eq. 19 approximately 
equals Up/us — m, so Sce is of order me. Ignoring higher orders of m for now, we can write 
= + 2e5eE, where the new e is independent of A' to order m. Now Eq. 14 becomes: 



15 iJ,'a 2 fi 15 jji'nsa 



16 r'-^ I 8 /xripr''' 
X [(1 + cos^ i) sin(2u;) sin(2A' - 20) + 2 cos i cos(2u;) cos(2A' - 20)] | 

x|(l + cos^i) cos(2cj) cos(2A' - 20) + 2 cosi sin(2c<;) sin(2A' - 20)|; (20) 

here the first {} arises from our changes in e. We seek the secular effects arising from evection, 
so we are interested solely in those parts of (20) that survive integration over A' from to 
27r; these are those terms that contain cos^(2c(;) cos^(2A' — 20) or sin^(2c<j) sin^(2A' — 20). 
Such terms, when averaged over A' (cf. Eq. 15), produce factors of |cos^(2c<;) or isin^(2c<;), 
respectively. Since we are interested only in the leading terms coming from evection, we can 
arguably ignore the eUipticity of the planet's orbit. Accordingly, the averaged i?^; is 

< ^(4)'(-)-aVyr^(l + cos2i)cosi[cos2(2a;) + sin2(2u;)]. 

The dependence on u obviously vanishes when the last two terms are added, and the lead 
factor [a/a''^Y can be given instead by the mean motions. Hence 

225 

< Re >= n^a^e^Vl - 6^(1 + cos^ i) cos i. (21) 



The extra uj (Danby 1992) arising from this term is 
Replacing UpUsa^/ii with m and grouping the terms, we get 



/dw\ ^ 225^^2n^ cos J[2(l - e^) - e2](l + cos^ i) + e^l + 3 cos^ ^)|. 
V dt / £/128//v J 



'du!\ 225 9 . n . 9, 

— = m n„cosz(2 - sm^z - e ). (22) 
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Similarly, 



(^] =-^mV'(4-3sm^i). (23) 
\ dt J E 128 ^ ^ ' ^ ' 



We note that Qe has as the leading term in eccentricity, whereas cj^; is non-zero even for 
circular orbits. In the case when i = and e = 0, we = + ^e = (225/32) m^np, which 
is the second term in Eq. 6, and comes from classical lunar theory. Therefore, it is clear 
that evection is by far the largest correction to any purely secular theory, such as Kozai's. 
Although Eq. 22 describes the only evection-related secular perturbation that affects circular 
orbits, it does not represent the full effect of the evection term in the disturbing function (Eq. 
14). To derive Eqs. 22 and 23 we used only the perturbation that a satellite's eccentricity 
suffers from evection. The fact that Eq. 14 depends on fl indicates that we should take into 
account that the inclination, too, suffers a periodic perturbation. 

The equivalent of Eq. 18 for inclination is: 

/di\ Ib^'usO? 1 

\dt)E~ 8 nr'^ Vl - sin i ^ 

x{(l + cos^ i)[cos{2uj) sin(2A' - 20) + cosisin(2a;) cos(2A' - 20)] 
-2 cos i [sin (2a;) cos(2A' - 20) + cosicos(2a;) sin(2A' - 20)]}. (24) 
Integrating as we did for Eq. 18, 

. 15 Us n'a^ 1 

dlE — o = X 

16 Up jj,r' V 1 — sm ? 

x{(l + cos^i)[-cos(2u;)cos(2A' - 20) + cosisin(2u;) sin(2A' - 20)]. 

-2 cos i[sin(2a;) sin(2A' - 20) - cos i cos(2a;) cos(2A' - 20)]}. (25) 

Since the inclination enters Eq. 14 only through cosi, it is useful to simplify the last expres- 
sion and convert it to ^(cosi)^: 

- , ns u'a^ 15 . o . 

olcostjE — T sm 

^ Up nr'^ 16 

x{cos(2a;) cos(2A' - 20) + cosisin(2a;) sin(2A' - 20)}. (26) 

Putting (26) back into (14) and integrating the resulting expression in the same manner as 
the term in Eq. 20, we get 

< Re' >= -— m nla^ ,^—= cos % sin^ % . (27) 

256 " ^ ^ 
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The resulting u and Q, caused by this term are 
225 



/dux 



sin^ i 



225 
128 



and 



/dQ\ 
\dt) E' 



l-e2 

pC0si(2e^sin^i — ^cos(2i)) 

225 o € ' 



(2 cos^ i — sin^ i)} 



256 



-m n. 



Pi 2 

1 — e"^ 



(2 - Ssin^i). 



(28) 



(29) 



Several characteristics of these terms are important. Both (28) and (29) are zero for circular, 
planar orbits, so they make no contribution to (6). Likewise, We' is zero for all orbits with 
i — 0, regardless of eccentricity Since the leading term in Eq. 28 is e^cosisin^i, this 
expression (for a satellite with e = 0.5 and i — 45°) is an order of magnitude smaller than 
the leading evection term (Eq. 22). 

The only term in Eq. 11 that wc have not addressed yet is Rj, or the "nodal evection" 
term. Since it does not depend on a;, it has no effect on the eccentricity, but only on the 
inclination (see the Lagrange equations in Danby 1992). The resulting perturbation is: 



diJi 



,3 f^^A^ - + T^'l «i"^sin(2y - 2Q). 



When integrated, this perturbation produces the following periodic term in sin^i: 



5(sin^ i)i — 



Ug jl'a^ 1 f3^^ 2\ 15 2i -2 



„ , ,-(1 - e^) H e^l cosisin^icos(2A' - 2f]). 



(30) 



(31) 



Once (31) is substituted into (13) and the resulting expression is integrated over A', we get 



< Rt >- 



9 1 9 

— mnia^—, (1 + 3e^ + -e^) cosi sin^ i. 
64 ^ Vl - 4 



(32) 



Using Lagrange's equations, we can calculate the secular effect of nodal evection: 



fduj 



9 



KdtJi 64 ^ 

1 + 



^ ^ 2 1 + 3e2 + (9/4)eS . . 2 . 
6 + 9e^ H ; \' ' cosz sm^ i 



(9/l)e^ 



l-e2 



cosi(2 cos^ i — sin^ 1) 



— 'm rip cos I 



9 2. 2 (l + 3eV2)2 ^ ■ 

3 + -e^ sin^ i - ^ cos 2i 

2 1 



and 



fdn^ 9 2 (l + 3eV2)% ^ 



(33) 
(34) 
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While the numerical coefficients in front of the nodal evection term are much smaller than 
that for apsidal evection, nodal evection still needs to be addressed if we want to determine 
the rate of nodal regression accurately. Eqs. 33 and 34 indicate that w is not affected by 
nodal evection if the orbit has no inclination, while there is always some contribution to 
VL. When e = and i = 0, fi/ = 9/32 m^np, which is identical to the second term in the 
expansion in m of the lunar nodal regression rate (the first being —3/4 mnp; Brouwer & 
Clemence 1961, see also Saha & Tremaine 1993). 

Until now, wc have determined the exact form of all terms of order m'^np or lower that 
enter the expression for the precessions of uj and Q. We have, however, neglected any terms 
containing higher powers of m, which are bound to arise in multiple places. Eq. 6 clearly 
demonstrates that w converges very slowly over powers of m, implying that many terms 
have to be included to achieve an acceptable accuracy for satellites with m's comparable 
to the Moon's (0.075). To illustrate this. Fig. 1 shows how the apsidal precession rate 
depends on the satellite's eccentricity in various theories and in direct numerical integration. 
The planet's eccentricity is taken to be very low (0.01) to avoid the effects arising from the 
octupole term (Section 4). Fig. 1 indicates that, while evection comprises about 2/3 of 
the discrepancy between the Kozai theory and the numerical integration, the accuracy of an 
analytical theory using only terms up to is limited, especially for low-eccentricity orbits. 

Even though the direct derivation of terms with higher powers of m appears to be very 
difficult if we are to use the same method as in this section, we can make some reasonable 
assumptions about their form. First, we note that all the terms in uj contained cosi, 
making them symmetric rather than antisymmetric with respect to the direction of orbital 
motion. Comparing our result with Eq. 6, we note that the variable m actually has a 
different meaning in our approach from that in lunar theory. In Eq. 6, m can be both 
positive and negative, to accomodate prograde and retrograde orbits, respectively. In the 
Kozai theory, just like in the present paper, the inclination is dealt with explicitly, and m (as 
well as the mean motion of the satellite) is always considered to be positive. So, in Eq. 6 the 
leading term changes sign for retrograde satellites, while the term is always positive. Eq. 
4 states that uj is always positive (at least for low-i orbits); however, owing to the different 
definition of zu for retrograde orbits {fl — uj), the motion of the line of apsides is negative. 
Likewise, the rri^ term in Eq. 6 is always positive, which requires the presence of a factor 
cosi in the cu term causing it (Eq. 22). We therefore propose that every subsequent term in 
UJ that contains m to a power n, should contain cos i to the power of n — 1, in order to make 
the terms alternately symmetric and antisymmetric. 

Figure 1 indicates that terms in w beyond the evection itself decrease with increasing 
e, so we need to address this dependence. We observe that the secular-plus-evection curve 
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intersects both the Kozai curve and the numerical values at approximately e = 0.8. Based 
on this, we will postulate that higher-order terms in m behave like (1 — 3e^/2). Note that 
-rij coming from evection has the same dependence on e in the planar case (Eqs. 22 and 
23). Unlike the dependence of higher-order terms on inclination discussed in the paragraph 
above, this approximation is more pragmatic and does not pretend to be exact. 

So our final equations for cu and Cl are 

/duj\ /duj\ /duj\ /duj\ /duj\ /duj\ 

\dt) = \li)K^\li)E^\ltt)E'^\dt)l^ \dt)L 

and 

fdVt\ /dVt\ /dVL\ /dVt\ /dVt\ 

where the u terms with subscripts K, E, E' and / are given by Eqs. 4, 22, 28 and 33, 
respectively. Likewise, Eqs. 5, 23, 29 and 34 define the Cl terms with respective subscripts 
K, E, E' and /. We define ojl as: 

, 10 

(^)l " ""^^^ ~ ''''' ^ ^^^^ 

where the Cj are taken from Tisserand (1894), and are very similar (but not identical) to the 
coefficients in Eq. Al of Saha & Tremaine (1993). To obtain a closed system of equations that 
can be evolved numerically, we also need to include expressions for the secular evolution of i 
and e from the Kozai theory (Eqs. 2 and 3). Evection induces only short-period variations 
in eccentricity which were fully accounted for during the averaging process, and has no 
consequences for the long-term evolution of the eccentricity (except for distant prograde 
orbits that are unstable; Nesvorny et al. 2003). The w predicted by Eqs. (35) and (36) for 
a planar, prograde case is plotted in Fig. 1, and its agreement with numerical integrations 
is quite satisfactory. In the next section, we compare the predictions of the present theory 
to the direct numerical integrations for a wide range of mean motions, eccentricities and 
inclinations. 



3. Comparison with Numerical Integrations 

Our first step in evaluating how closely our analytical model represents the behavior 
of real satellites is to compare it to a set of directly integrated test particles. In order to 
classify possible deviations of the model from the observed dynamics, we will first look at 
orbits that either have high eccentricity or high inclination, but never both. Only after that 
will we apply the model to the real sateUites, which often exhibit both. Recall that in the 
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conventional Kozai theory, (1 — e^)^/^ cosi (the orbit's angular momentum that hes normal 
to the reference plane) is conserved (see Holman et al. 1997, Carruba et al. 2002). Thus 
changes in e and i are coupled. Hence, throughout this chapter, we will characterize bodies 
by their "minimum inclination" and "maximum eccentricity". By this we are referring to 
the values reached at one extreme of the Kozai cycle, when u = 90°. This way the orbits 
are uniquely defined in a secular model which ignores octupole and higher terms, since it 
depends only on a, e, i and u (the inclusion of the octupole term brings in a dependence on 
w; Yokoyama et al. 2003). We prefer these "extreme elements" over the mean ones since 
two irregular satellites with quite different orbits can have the same mean elements. Also, 
these extreme elements are used as initial conditions in all the integrations shown in this 
chapter. If the Kozai resonance is possible for a particle's combination of a, e and i, an 
initial uj = 90° will put it in the libration region. This way, all of our test particles above 
the critical inclination will have librating uu, enabling us to find the exact location of the 
boundary. 

The first group of integrations deals with the secular behavior of test particles having 
very low minimum inclination (5°, or 175° for retrograde cases) and maximum eccentricities 
varying from to almost 1. In all secular theories, Eqs. 2 and 3 were used for the evolution 
of inclination and eccentricity, respectively. All of the continuous lines in the figures were 
obtained through advancing the relevant secular equations by a Burlisch-Stoer-type numer- 
ical integrator over multiple precession periods, and then computing the average precession 
rate. The discrete points show the results of a direct numerical simulation for the same 
parameters. Just as in Fig. 1, the eccentricity of the planet ( "pseudo- Jupiter" ) in the sym- 
plectic intergration was taken to be only 0.01 to avoid any interference from the octupole 
term (shown in Sec. 4 to be proportional to Cp) . 

It is clear from Fig. 2a that our model is a very good approximation to the secular 
behavior of a prograde, low-i satellite for all values of eccentricity. While the addition of 
higher-order terms (Eq. 37) was needed to accurately describe the motion of the line of 
apsides, evection alone is sufficient for a satisfactory description of nodal regresion for all 
e. Fig. 2b shows that our model of apsidal precession gives excellent results for retrograde 
satellites, too. However, the precession rate of the node differs from the one predicted for 
particles with e > 0.5 (although our model is still more accurate than Kozai's as long as 
e < 0.9). The most hkely cause of this disagreement is the existence of higher-order terms 
in m that affect the nodal precession. This hypothesis is supported by Fig. 2c, in which the 
theory and numerical simulation of the nodal precession noticeably diverge already at e > 0.4. 
In contrast, our predictions for the apsidal precession are still reasonably accurate, and the 
slight shift between two curves is most likely a consequence of Cj's with j > 10, that are 
not included in expression (37). Unfortunately, since these "phantom terms" in fl vanish if 
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either e = or i = (but not when i — 180°), there is no easy way to synthesize them on the 
basis of lunar theory, as we did with apsidal terms in Eq. 37. This inabihty of our empirical 
approximation to describe accurately the nodal motion of distant and eccentric retrograde 
orbits is irrelevant for the discussion of secular resonances involving only the longitude of 
pericenter (Sections 4 and 5). On the other hand, as we shall see, the calculation of the 
boundary of Kozai resonance is affected, with the model giving us erroneous results for the 
secular behavior of at least one known irregular satellite (see our later discussion of Fig. 4b). 

Figures 3a-c test the behavior of our model with changing minimum inclination. The 
maximum eccentricity in all three panels is 0.2. In all three panels wc sec that at a certain 
inclination the average precession rates of the line of apsides and of the line of nodes converge 
and stay equal to each other at all higher inclinations. When those two rates are identical, 
the argument of pericenter hbrates around 90° or 270°; this behavior is usually known as 
the Kozai resonance (see Sec. 1). It is interesting that the Kozai resonance is not achieved 
by both secular solutions at the same inclination. In Fig. 3a the Kozai resonance is reached 
by our model at a minimum inclination of about 47°, while the Kozai model, as expected, 
predicts the boundary to be at 39.2° (in Kozai theory this result is independent of the mean 
motion). The numerical simulation, however, agrees closely with our model (the observed 
slight divergence within the libration region arises because the points actually show the 
average of inclinations at which uu passes through the libration center, regardless of the 
direction). This higher threshold for the Kozai resonance among distant prograde objects 
led to erroneous early predictions based on purely secular models that Siarnaq (formerly 
S/2000 S3) should be in Kozai resonance (Vashov'yak 2001); direct numerical integrations 
showed it to be in the secular resonance instead (Cuk et al. 2002a; Nesvorny et al. 2003). In 
Figs. 3b and 3c, the Kozai resonance is reached at inclinations larger (i.e., closer to 180°) 
than 140.8° by both the numerical simulation and our model. Although our model puts 
this threshold at slightly larger inclination than the simulation, it still approximates the 
numerical results much better than does the traditional Kozai model. 

Figures 2 and 3 suggest that our model fairly accurately predicts the precessional rate 
of the line of apsides for all the cases tested, and that it can predict the motion of the line of 
nodes for all prograde orbits and for those retrograde orbits that have low eccentricity. Figs. 
3a-c show that the boundary between a;'s libration and circulation is not fixed but varies 
with different mean motions. The shift in this boundary reduces the range of inclinations 
over which prograde orbits can be in Kozai libration, while it expands the range of retrograde 
inclinations that allow libration. 

To identify the region in which Kozai resonance is possible, we determine the inclination 
at which circular orbits switch from circulation to libration as a function of the orbital period 
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(measured by m). To make the plot continuous, we define inclination to be always smaller 
than 90°, with m taken to be positive for prograde bodies and negative for retrograde ones. 
The boundary is calculated by substituting e = and cu — 90° in the separate terms in Eq. 
35, which is then equated to zero: 

dw ^ r . 2 • 225 ., . 2 A 
— — 2 — 5 sm I + -^-j-^ cos z(2 — sm t) 

9 ^ _ _ 

H m cos i(3 sin^ i — cos(2i)) + Cj ^ cos^ ^ i = . (38) 

32 . 

3=3 

The classical result that Kozai resonance occurs for i > 39.2° (or i < 140.8°) comes from the 
first two terms on the right (sin^i = 2/5). Equation 38 is solved iteratively for i, by first 
assuming that sin^ i — 0.4, and then finding successive approximations to i using 

. r2 45 . . 2 . X 

in+i — arcsmj- + —m cos in[2 — sm z„j 



9 1 ^° _ 

+— mcosi„[3sin^i„ - cos(2i„)] + - ^cos-' ^inV^'^- (39) 

i=3 

This expression converges well for m values that are typical of the known irregular satellites 
(-0.2 — ^ 0.1); since m and cosi always appear together in Eq. 39, the expression does not 
need to be modified now that we switch to a new convention in which retrograde orbits have 
m < and i < 90°. Fig. 4a plots the location of the Kozai resonance as defined by this 
procedure. 

Since circular orbits, strictly speaking, should exhibit no Kozai cycle in inclination over 
the precession (or libration) period of uu, this inclination is both a minimum and an average 
at the same time. So we also plot, as individual points, the m's and minimum inclinations of 
some known irregular satellites that lie in the same region. These minimum inclinations were 
determined on the basis of a 300,000-year symplectic integration for each sateUite. For these 
integrations we used our home-made orbital integration software, which follows the stan- 
dard algorithm devised by Wisdom & Holman (1991). The initial conditions were obtained 
from JPL's Horizons ephemeris service on December 29, 2003, in the form of the osculating 
elements for the four giant-planet barycenters and the irregular satellites themselves for mid- 
night, September 24, 2003. The minimum inclination is defined as the average inclination to 
the planet's orbital plane for times when coso; < 0.1 and the eccentricity is larger than the 
average value for the integration. Different symbols represent the sateUites of each of the 
four giant planets (see figure caption); those sateUites that are found to have hbrating a;'s 
are also marked by large triangles, as well as identified. The satellites of Jupiter that were 
discovered in 2003 are not included in this plot, as their orbits were not well constrained 
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when the initial conditions were generated, although they are very likely to fall into the al- 
ready known orbital clusters (Sheppard & Jewitt 2003). An exception was made for S/2003 
J20, which is very unlikely to belong to any of the known groups, and is therefore probably 
relevant to the question of the overall orbital distribution of the originally captured bodies. 

Two features of Fig. 4a require comment. First, clearly our theoretical prediction for 
the circular case poorly describes the behavior of the real satellites. More Kozai librators 
are foimd slightly below the line than above it; also, a couple of Neptunian retrograde 
irregulars that arc above the line are not thought to be Kozai librators. The reasons for 
this discrepancy are that almost all of the librating satellites have substantial eccentricities, 
and that the boundary between the two types of secular behavior shifts with increasing 
eccentricity. Just from the distribution of known librators it is possible to speculate that the 
boundary moves down for more eccentric prograde moons, and up for eccentric retrogrades. 
The only satellite for which this low-e approximation is accurate is Euporie (formerly S/2001 
JIO), which indeed docs have low eccentricity {cmax = 0.17). The other important feature 
of this plot is that the boundary goes through a local minimum in inclination close to 
Euporie's position (m = —0.13); for more distant retrograde orbits, the circulation-libration 
transition happens at increasingly larger inclinations. This change in behavior apparently 
occurs because, at more negative m's, terms of order 3 and higher in Eq. 39 become more 
important than the pair coming from evection. 

To account fully for the effects of eccentricity on the position of the Kozai boundary, we 
need to carry out a full integration of the equations affecting cj's motion (Eqs. 2, 3 and 35). 
Although Eq. 36 is decoupled from the other three {Q cannot influence other elements in 
the quadrupole-limited theory), we evolve it alongside the other three for comparison. We 
do this for three sets of test particles that form a dense grid in m and minimum inclination 
(the cell size is 0.01 xl°), with the sets having maximum eccentricities of 0.2, 0.4 and 0.6. 
The planet's parameters are those chosen for Figs. 2 and 3. After the particles' orbits are 
advanced over tens of precession periods with a Burlisch-Stoer-type integrator, the average 
precession periods of uj and Q are computed for each particle. Then, for each m, we seek 
the particle with the lowest inclination for which uj < O.IQ (a relative criterion is required 
since close-in moons always have longer secular periods than the more distant ones), and we 
use this to indicate the onset of the Kozai resonance. These lowest minimum-inclinations 
for which cu is librating are connected by continuous lines in Fig. 4b. 

As expected, the boundary for the Kozai resonance among the prograde satellites moves 
to lower inclinations with increasing eccentricity. Now all of the known prograde librators 
are either above, or less than a degree below, the boundary for their eccentricity range (the 
way that boundaries were computed can lead to critical inclinations being overestimated by 
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as much as one degree). No circulating satellites are found to have minimum inclinations 
above the boundary. We conclude that our model predicts the position of the separatrix 
between circulation and libration among prograde moons accurately enough for the purpose 
of discussing their overall distribution and behavior. It is also noteworthy that all but one 
(S/2003 U3, the orbit of which still has to be confirmed) among the prograde librators are 
fairly close to the circulation-libration boundary, which may offer some clues about their 
origin and past evolution (Section 6). 

On the other hand, the agreement between the Kozai boundary predicted by our theory 
and the behavior of high-i retrograde satellites is poor, especially for the more distant bodies. 
Instead of shifting to higher i as e increases (as might be suggested by the positions of the 
real satellites on the plot), the computed boundary appears to do exactly the opposite. The 
cause of this discrepancy is clearly illustrated in Fig. 2c. Our model fails to predict the nodal 
precession rate of the distant and eccentric retrogrades with any accuracy, while it does much 
better when applied to the motion of the line of apsides. Since the Kozai resonance deals with 
the behavior of the argument of pericenter (defined as the angle between the ascending node 
and the periapse), our model should be expected to produce major errors when trying to 
predict its location. Fig. 2c indicates that, according to our model, the precession of the line 
of nodes (of a distant, retrograde satellite) slows down with increasing eccentricity, while 
the numerical simulation suggests that it actually accelerates. Since the Kozai resonance 
happens because of a reduced apsidal precession rate with increasing inclination (Fig. 3c) 
and, ultimately, its synchronization with the nodal precession rate, the inclination at which 
this synchronization arises must depend on the speed at which the nodes move. Fig. 3c 
suggests that a faster nodal precession would lead to the Kozai resonance starting at a 
higher inclination (by "higher" here we mean "closer to 90°"), and vice versa. Therefore, 
the lines plotting the circulation-libration boundary for retrograde satellites with higher 
eccentricities should be consistently above the zcro-e curve shown in Fig. 4a. So the upturn 
seen at m = —0.13 in Fig. 4a should happen at higher i and be even more pronounced for 
larger e. It is easy to anticipate that for a maximum eccentricity of about 0.7, the boundary 
should pass between the positions of the new Neptunian moons S/2002 N4 (labeled on the 
plot) and S/2003 Nl (the box directly below N4), which have similar orbits except that the 
former is thought to be a librator whereas the latter's u circulates. Of course, the present- 
day orbits for these objects might differ from the true ones, but this does not make the 
numerical simulation which predicts that these two sets of initial conditions lie on opposite 
sides of the boundary be any less accurate. Having this in mind, we remark that the only 
two retrograde moons that are thought to be Kozai librators are also not far from the true (if 
not the modeled one) circulation-libration boundary, just like their more numerous prograde 
cousins. 
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We now identify the position of the secular resonance, defined as 'dj — 'ti^pianet — 0. This is 
shown in Fig. 5, which is generated from the same set of integrations as Fig. 4b. The shape 
of the secular resonance was found by locating, for each set of particles with the same m and 
e, two test particles on neighboring inclinations that have their lines of apsides precessing in 
opposite directions (where w = Q + u for prograde, and w = Q — uj for retrograde particles). 
The critical inclination for the secular resonance was then found as an intermediate value 
between the inclinations of this pair at which w — (the dependence of on i was assumed 
to be locally linear). Strictly speaking, the condition that ~ is a necessary, but not 
sufficient, for the secular resonance's existence, and the apsidal locking itself cannot be 
described by a purely quadrupole theory. So the lines in Fig. 5 should only be considered 
as potential locations of the secular resonance, rather than the places where apsidal locking 
always happens. The feasibility and strength of this apsidal locking is discussed in the next 
section. Here wc will merely note that the positions of the four resonant (or pseudo-resonant) 
objects - Pasiphac {cmax = 0.50), Sinope {cmax = 0.32), Stephano {cmax = 0.32) and Siarnaq 
i^max = 0.49) - are in good agreement with the location of the secular resonance computed 
from our model. Small differences between our theory and direct integration are expected 
since we have ignored the octupole terms that are vital for the dynamics of resonant objects 
(Yokoyama et al. 2003). A possible correlation between the secular resonance's position and 
the overall distribution of the known irregular satellite groups is discussed in Section 6. 



4. Dynamics of Secular Resonance 

Figure 5 locates those positions in orbital-element space where the precession rate of the 
longitude of pericenter is zero, as computed by a numerical integration of Eqs. 35 and 36. For 
some time it has been known that the apsides of a slowly precessing irregular satellite's orbit 
can become locked in a resonance with the planet's perihelion (Whipple & Shelus 1993; Saha 
& Tremaine 1993). Figure 6 displays the evolution of the resonant argument \l/ = w — Wpianet 
for a) Pasiphae, b) Sinope, c) Siarnaq and d) Stephano. These 300,000-yr simulations were 
carried out using a symplectic integrator, with the same initial conditions as in the previous 
section. It is obvious that the medium-term behavior of the resonant argument is different 
for each of the four objects. Pasiphae's resonant argument librates around ^ = tt for the 
integration's entire length. Sinope's pericenter alternates between circulation and libration; 
the center of libration is the same as Pasiphae's, but the amplitude is much larger. Siarnaq 
also exhibits an elaborate but periodic mixture of circulation and large-amplitude libration, 
only now the center of libration is ^ = 0. Finally, Stephano's pericenter simply circulates, 
albeit with a very long period (nearly 200,000 yrs). 
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This diversity of behavior, including the hbrations, cannot be explained by a quadrupole- 
only theory (such as ours), since the quadrupole term has no dependence on the orientation 
of the planet's line of apsides (cf. Eq. 16). Yokoyama et al. (2003) have derived the next 
two terms (in a/ a') of the disturbing function for satellites with arbitrary e and i, and have 
found that the octupole term (i.e., the one that contains {a/a'Y) is responsible for apsidal 
locking. The complete octupole term (Yokoyama et al. 2003) is: 

Ry = 7^ [(—3 + 33 cosi + 15 cos^ i — 45 cos^ i)bi cos(zusun — + 
64 

+(—3 — 33 cos i + 15 cos^ i + 45 cos^ i)bi cos{zusun — ^ — uj) 
+(15 — 15 cosi — 15 cos^ i + 15 cos^ i)b2 cos(ro5„„ — Q + 3a;) 
+(15 + 15 cos i — 15 cos^ i — 15 cos^ i)b2 cos(ti75„„ — Q — 3a;)] , (40) 



where 



' 3 

= ^jrzr^^ (41) 



p 

.3 



and 



fei = -(5/2)e- (15/8)e^ (42) 



-(35/8)e^ (43) 



Yokoyama et al. were the first to be able to derive this octupole term for arbitrary e and 
i, but a simple version of this term has been featured in lunar theories. If we substitute 
cosi = 1 (i.e., i = 0) in Eq. 40 and ignore powers of e and Cp beyond linear, all the terms in 
Eq. 40 vanish except for the second, which becomes 

Ry,m = -^^^eepcosiwsun - ^ - u;) ^ -^nJa^^eepCOs(ro - wsun)- (44) 

Brouwer & Clemence (1961) mention expression 44 as the most important term in the lunar 
disturbing function that produces a perturbation independent of m. This is to be expected, 
since we have seen in Sec. 2 that m arises when short-period perturbations are averaged, 
while the octupole term is purely secular. Apart from this interesting feature, the octupole 
term is a very minor perturbation to the lunar motion, because the factor a /a' is very small 
(2.5 X 10-^) for the Moon's orbit. 

Yokoyama et al. (2003) correctly point out that the octupole terms become important 
only if the argument of one of their cosine functions changes very slowly. Near secular 
resonance, the arguments in one or the other of the first two terms actually do become 
very slow for retrograde and prograde objects, respectively (remember that zu — + cu ior 
prograde orbits, but zu — — cu ior retrograde ones). In order for the arguments of the 
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third and fourth terms to become stationary, a sateUite must have a very slowly precessing 
argument of pcriccnter, as well as a sluggish Such a resonance may actually be possible 
close to the boundary of Kozai libration, but these terms are irrelevant for our discussion 
of the secular resonance (for a detailed analysis of the dynamics close to the circulation- 
libration boundary, see Carruba et al. 2004). Therefore the portion of the octupole term 
relevant for the secular resonance (Yokoyama et al. 2003) is 

Rr = -^(-3 - 33| cosi| + IScos^i + 45| cos^i|) 6iCos*, (45) 

where it is implied that ^ = i7 + — zup for progradc, and ^ = Q — cl^ — zup for retrograde 
objects (where Wp = Wsun + vr is the planet's longitude of pcriccnter). The eccentricity and 
inclination of the satellite that appear in Eq. 45 represent average values over many secular 
(i.e., Kozai) periods, since zi; << ci;, by definition, for objects in the secular resonance. 

We will now describe how Eq. 45 may explain some of the diverse behavior among the 
known resonant and near-resonant satellites. First, we point out that the relevant parameter 
for resonant-locking is not the absolute magnitude of the octupole term, but how it compares 
to other perturbations, most notably the quadrupole term. Eqs. 12-13 show that all parts 
of the quadrupole term contain the tidal factor /x'a^/a'^, which effectively determines its 
strength. So, to assess the resonant term's relative importance, we divide Eq. 45 by ji'a^/a'^. 
Assuming Sp « 1, the resulting expression is 

K = -ep^S {cos i)hi (e) cos , (46) 
a 

where S{cosi) is 

S{cosi) = (-3-33|cosi| + 15cos2i + 45|cos^i|)/64. (47) 

We see that the resonant term is proportional to the planet's eccentricity, meaning that it 
will vary in strength as tp changes. Fig. 6b suggests that Sinope's episodes of libration and 
circulation alternate with a period of roughly 50,000 yrs. We recall that secular changes in 
Jupiter's eccentricity have a similar period (Murray & Dermott 1999). We therefore suggest 
that resonant-locking is drivern by Cp. To illustrate this correlation between the behavior 
of Sinope's resonant argument and Jupiter's eccentricity. Fig. 7a gives Jupiter's Cp versus 
Sinope's ^ during the course of the numerical integration shown in Fig. 6b. It is clear that 
the region close to ^' = is "forbidden" (and therefore circulation is impossible) during those 
times when tp > 0.045. While the overall character of Siarnaq's ^ appears very different 
from Sinope's, Fig. 7b shows that there is a correlation between Siarnaq's and Saturn's 
Cp, too. For the duration of our integration, Siarnaq's ^ follows a periodic trajectory which 
passes through ^ — n only when Saturn's eccentricity is at its minimum. This apparently 
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happens during every third minimum in Cp, or about every 150,000 yrs. At all other times 
Siarnaq's pericenter avoids anti-alignment with Saturn's. So we conclude that the medium- 
term changes in the behavior of a resonant satellite's ^ appear to be forced solely by the 
secular variations in the eccentricities of Jupiter and Saturn. A longer-period variation has 
also been observed (Nesvorny et al. 2003) but its causes are likely to be much more complex 
(see below). 

We can use Eq. 45 to investigate the stability of apsidal alignment versus anti-alignment 
for different orbits. Since &i(e) is always negative (Eq. 42), the sign of R'^ depends only on 
S'(cosi) and cos ^. To determine which point will be a center of oscillations in tu, we need to 
obtain an equation of motion for w near the resonance. In secular resonance, zu is dominated 
by the quadrupole term, while the long-term e is caused solely by the octupole term (Lee & 
Peale 2003). Since the precession of w caused by the quadrupole term has a zero close to 
the exact resonance, we can linearize zu{e) around that point: 

ZD = K^{e- Cres), 

where e^es and are constants. By differentiating this expression with respect to time, we 
get: 

(PvD de 

which, through Lagrange's equations (Danby 1992), becomes 

d'^uj ^ nay/l — cosi dR 
dt"^ ^ //e I cosi| dw 

The sign function (i.e., cosi/| cosi|)is needed since = Q -|- a; for prograde bodies, but 
= Q — a; for retrograde ones. The only term in the secular disturbing function that 
depends on w is the octupole, so we can write: 



(^}w cosi na\/l — cy ^, .\ , . ^ 

— - = — -K^ — S{cost) 6ism*. 

dt'^ I cos z I lie 64 

If we consider that Wp is constant, this can be re-written as: 

d^^ cos i 

-n^ - Qp-, ttK^ S{cos i) sin ^ = 0, (48) 

|cosz| 

where Qp is a positive quantity, independent of zu. Eq. 48 is the well known pendulum 
equation (Murray & Dermott 1999). Depending on the direction of orbital motion and the 
signs of and 5'(cosi), the stable center of ^ oscillations can be either or tt (in the case 
of negative and positive product of those three terms, respectively). While is a complex 
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function of orbital elements and must be obtained through numerical integrations (cf. Saha 
& Tremaine 1993), S'(cosi) depends only on inclination. Fig. 8 plots how S'(cosi) varies 
with average inclination. S is positive for i < 41°, and negative for higher inclinations. The 
positions of the four (near-)resonant satellites are marked: Pasiphae and Sinope (as well 
as Stephano) are in the region where apsidal anti-alignment for retrograde bodies requires 
< 0, whereas prograde Siarnaq's stable point of ^' = requires > for its set of 
orbital elements. In the case of Sinope, our numerical simulations show that the episodes of 
negative circulation of ^ correlate with higher e, which indeed implies that < (see Fig. 
14 and its discussion in Section 6). Similarly, Fig. 17 of Nesvorny et al. (2003) shows how 
negative circulation for Siarnaq coincides with smaller eccentricity, so > for that moon, 
as expected from Eq. 48. Therefore, Eq. 48 clearly explains why the libration center for 
Pasiphae and Sinope is diametrically opposite from Siarnaq's. Although the transition from 
one libration center to the other is related to the direction of the satellite's motion (through 
the different definitions of zo for prograde and retrograde objects), retrograde satellites could, 
in principle, exist with orbits librating around ^' = 0. However, Fig 5 suggests that secular 
resonance is possible only for prograde bodies with inclinations above ~ 40°, and for retro- 
grade satellites with inclinations below that value, making the product cosi S{cosi) almost 
always negative. Barring some unexpected behavior of K^{a, e,i), any new retrograde res- 
onators that are found are likely to be in the "Pasiphae- regime" (librations around ^ = tt), 
while the prograde ones can be expected to be in the "Siarnaq-regime" (with \1' = as a 
stable point). Exceptions would be resonators with inclinations of about 40°, at which all 
three relevant factors in Eq. 48 change sign, but resonant locking is unlikely to be strong for 
them in the first place, due to their smaller a/a' ratio (see below). 

In order to probe the cause for the very different behaviors among the four objects 
mentioned above, we will need to take into account all the terms in Eq. 46. Table 1 lists 
a/ a', the average e, &i(e), the average i, S{cos i) and the relative strength of the resonant term 
for each of the four resonant moons; averages are taken over the length of the integration. 
The last quantity was obtained by dividing R'^ for the each sateUite by Pasiphae's E!^. Here 
we assumed that the instantaneous eccentricities of Jupiter, Saturn and Uranus are equal to 
each other (they all vary between and 0.1 over secular timescales). We see that Sinope's 
resonant term is of the same order of magnitude as that of Pasiphae (133% of Pasiphae's), 
while those for Siarnaq and Stephano are one and two orders of magnitude smaller (19% 
and 2%, respectively). This comparison clearly shows why Stephano's pericenter cannot 
currently get locked into a resonance with that of Uranus, as its resonant term is much too 
weak. The low value for Stephano's resonant term mostly results from its small a/ a' ratio, 
although its eccentricity and inclination are also less conducive to resonance than is the case 
for the other three moons discussed here. In turn, the small a /a' ratio comes not only from 
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the relatively close orbit of Stephano (m=0.022), but also from the small size of Uranus's 
Hill sphere. The size of the Hill sphere (in terms of a') depends on the planet's mass as 
^/i/3^ making the Hill spheres of Uranus and Neptune significantly smaller fractions of their 
a' compared to those of Jupiter and Saturn, m = [a/ RuY^"^ measures the size of a satellite's 
orbit relative to the planet's Hill sphere, so two moons that have the same m but that orbit 
different planets can have very different (a/a')'s (for example, Himalia and the Moon have 
comparable m's, but Himalia's a /a' is an order of magnitude larger). 

But, if Sinope's resonant term is larger than Paisphae's, why is Sinopc only occasion- 
ally librating, while Pasiphac appears to be deep within the libration region? The reason 
must lie in the initial conditions, which determine not only the libration"s amplitude, but 
if librations are possible in the first place. It is tempting to conclude that Sinope's larger 
libration amplitude tells us something about the origin and evolution of these two moons. 
However, we must remember that integrations shown in Fig. 6 cover only 300,000 yrs, an 
insignificant fraction of the Solar System's lifetime. In order to make any inferences about 
the intrinsic differences between these moons' orbits, we need to examine their behavior 
over much longer times. Nesvorny et al. (2003) have conducted numerical simulations of 
the known irregular satellites' orbits over 10^ yrs. They conclude that both Pasiphae and 
Sinope show intermittent resonant and near-resonant behaviors, with the switching between 
the two regimes occuring on 10^-yr timescales. For example, their Fig. 17 suggests that, in 
4 X 10^ yrs, Sinope will be a strict hbrator while Pasiphae's ^ will circulate. Judging from 
these results, it is very likely that the resonant arguments of both Pasiphae and Sinope are 
chaotic on timescales shorter than the age of the Solar System, precluding us from obtaining 
any direct clues as to their origins from their present libration amplitudes. 

Nesvorny et al. (2003) also found that Siarnaq occasionally exhibits relatively short 
episodes of pure libration. These episodes never last longer than 5 x 10^ yrs, unlike the 
more stable states of Pasiphae and Sinope. This result is consistent with our estimate that 
Siarnaq's resonance is about an order of magnitude weaker than those for the two Jovian 
resonators. Additionally, Nesvorny et al. do not report any resonant behavior for Stephano's 
which is in line with our conclusion that Stephano's resonant term is very weak (Table 1). 

Based on these results, as well as those of Yokoyama et al. (2003) and Nesvorny et al. 
(2003), we can make some general conclusions about the viability of the secular resonance 
for different orbits. Resonant-locking is very unlikely for the orbits in the middle portion of 
the continuous lines in Fig. 5, as the (a/a')'s of such orbits are small and their inclinations 
put them close to the zero of S{cosi). Librations should be more common for the resonant 
satellites along the edges of Fig. 5, especially among the distant retrograde moons. Also, 
librating behavior is more likely to be seen among the Jovians and Saturnians than the 
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Uranians, due to the relatively smaller Hill sphere of Uranus. Any present-day ^ hbrations 
among the Neptunians are highly unlikely, given the very low present eccentricity of Neptune. 
A satellite with a confirmed near-resonant orbit around Neptune (S/2002 Nl is the best 
candidate so far) might be a remnant of past locking, and therefore may indicate a higher 
primordial eccentricity of Neptune (see Sec. 6). 



Despite these benefits of secular models, they can never describe all the phenomena 
we see among the irregular satellites. For example, since they are averaged over the mean 
motions of the satellite and the planets, they are unable to register resonances that involve 
mean motions, which are known to be important for irregular satellite dynamics (Saha & 
Tremaine 1993; Nesvorny et al. 2003). We here report the serendipitous discovery of one 
such resonance, that involves the Great Inequality of Jupiter and Saturn, and that can have 
a suprisingly strong effect on the orbits of Saturnian irregular satellites. Recently, Carruba 
et al. (2004) have independently found several related resonances for Saturnian irregulars, 
although involving different resonant arguments from that discussed in this section. 

Fig. 9a shows the evolution of eccentricity for a test particle orbiting Saturn. The 
initial conditions for the particle are a = 0.1 AU, e = 0.4, i = 37.5°, u = 90°, Q = 180° 
and M — 45°, while the four giant planets start the integration with conditions identical to 
those at midnight, September 24, 2003 (see Sec. 3). This orbit should lie deep within a zone 
where satellite orbits are stable (Nesvorny et al. 2003). However, as Fig. 9a clearly shows, 
the eccentricity of this test particle exhibits large and irregular variations, with the average e 
growing from 0.3 to almost 0.5 in the course of our integration. Given that the time span of 
300,000 yrs is just a small fraction of the total age of the Solar System, we have to conclude 
that such an object would probably be unstable on longer timescales, with a hkely fate of 
colliding with lapetus (such a coUision becomes possible once e > 0.75). Integrations of a 
larger set of bodies indicate that the relevant parameter for this instability is the precession 
period of the longitude of pericenter. Fig. 10 plots a measure rj (defined below) of the 
eccentricity variation against the w rate for 440 test particles near the one shown in Fig. 
9a. These integrations spanned 30,000 yrs each and the range of initial conditions was: 
a = 0.1 - 0.13 AU, e = 0.2 - 0.875, i = 30° - 45°, with the steps being 0.01 AU, 0.075 and 
1.5°, respectively (for all bodies uj = 90°). The parameter rj is defined as 



5. Effects of the Great Inequality 
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where 
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and 




edt. 



A high value of r] for a test particle indicates that e has either a long-period (> 3000 yrs) 
oscillation, or a secular trend. Fig. 10 clearly shows the close correlation between the 
particles' rj and uj. The high-?7 particles tend to have ~ or ~ 5.4 x 10~^. The first of 
these two features can safely be identified with the secular resonance (Sec. 4). The apsidal 
precession period associated with the second feature is about 1860 yrs, which is close to twice 
the period of the Great Inequality of Jupiter and Saturn (883 yrs). The Great Inequality 
is a consequence of the mean motions of Jupiter and Saturn being close to, but not in, 5:2 
commcnsurability. This near-resonant perturbation has a significant effect on the orbits of 
both planets, and it prevented many early reascarchcrs from completely harmonizing theories 
based on Newtonian gravity with the observed planetary motions (Baum & Sheehan 1997). 
The 883- year period corresponds to the residual motion of the near-resonant argument 5Xs — 
2\j. In order to make a valid term in the disturbing function, consistent with d'Alembert's 
rule (which requires the sum of numerical coefficients to vanish), this argument has to be 
complemented by at least three other angles, with negative coefficients. The longitudes of 
the lines of apsides and nodes of Jupiter and Saturn (which change slowly) commonly fill this 
role, so a number of observed strong perturbations of these planets' orbits have periods close 
to 900 yrs. Based on all this, it is straightforward to conclude that if we subtract twice the 
zu frequency of our high-?] test particles, the resulting argument would change exceptionally 
slowly. Since half of the apsidal precession period of those test particles is still longer than 
the 883-year Great Inequahty period, we conclude that the third additional term should be 
a prograde secular angle, with vos being the obvious choice. So we define our candidate 
resonant angle ^ as 



Fig. 9b plots the evolution of ^ for the test particle whose eccentricity history was shown in 
Fig. 9a. Not only does ^ hbrate at times, but the correlation between the behaviors of e and 
^ is impressive. Apparently, the largest change in the eccentricity occurs when ^ librates, 
whereas e has much less of a secular trend at times when ^ circulates. On the time scale of 
this integration, the evolution of ^ appears to be chaotic, although there yet may be some 
longer-term periodicities. In any case, we feel that the identification of this perturbation 
with the Great Inequality, and the specific resonant argument ^ is firm. We should also note 
that Fig. 10 suggests that no other harmonics of the Great Inequality appear to affect the 



^ = 5A5 — 2Aj — 2w — ZDs- 
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test particles appreciably, at least within the range of apsidal precession frequencies found 
among our set of test particles (cf. Carruba et al. 2004). 

We now address the circumstances (combination of orbital elements) for which the Great 
Inequality resonance occurs. Formally, it should happen for any bodies orbiting Jupiter and 
Saturn whose w circulation periods are around 1800 yrs. We found that, while detectable, 
this resonance does not strongly perturb Jovian satellites (even if wj is substituted in place 
of Ws in Eq. 50), or those retrograde Saturnians for which relation (50) holds (this is possible 
for some orbits with inclinations slightly above the secular resonance). Effects as large as 
those displayed in Fig. 9a arc restricted to the progradc Saturnians. Saturn's satellites are 
likely favored due to the larger effects of the Great Inequality on that planet's orbit, but the 
reason for the prograde bias is likely to be more complex. Fig. 11 locates the Great-Inequality 
resonance in m — imin space. The top three continuous lines plot the resonant location for 
different eccentricities (see caption). It is not surprising that no satellites at the present 
time inhabit the resonant region (the proximity of Siarnaq to the top line is irrelevant, as 
that line applies to bodies of much lower eccentricity). It is even more interesting to explore 
how the location of the resonance shifts with the changing period of the Great Inequality. 
Since the latter measures the distance of Jupiter and Saturn from the exact 5:2 resonance, 
even relatively small changes in the orbital periods of those two planets can lead to a large 
variation in the Great-Inequality's period. Since Saturn probably migrated outward and 
Jupiter inward in the early Solar System (Hahn & Malhotra 1999), they should have been 
further from the resonance in the past (since hng — 'irij > 0). The bottom solid line in 
Fig. 11 plots the location of the Great Inequahty resonance for the hypothetical epoch when 
(5A5 — 2Aj)~^ = 500 yrs (we chose Cmax = 0.6, which is a good approximation for most 
prograde Saturnian irregulars). With this change of parameters, the resonance moves to 
signficantly lower inclinations, approaching the elements of Albiorix, which presently has 
~ 820 yrs; so a somewhat shorter period of Great Inequality (~ 400 yrs) would definitely 
affect it. Such a shift in the Great-Inequality's period requires a change in Saturn's period 
of less than 1%, and Aas < 0.1 AU. Migration on that, or larger, scale is hkely for a wide 
range of initial parameters of the primordial planetesimal disk (Hahn & Malhotra 1999). 

Since Albiorix is the largest member of a satellite cluster that also includes Erriapo 
and Tarvos (Gladman et al. 2001; we will call it the "Gauhsh cluster"), the history of 
this resonance may be coupled to the origin of the whole group. It is unclear whether this 
is a coUisional group (Nesvorny et al. 2003) and, if it is, if the passage through the Great- 
Inequahty resonance happened before or after the coUisional break-up. Nesvorny et al. (2003) 
find that the velocity dispersion of the cluster is too large (30-60 m/s) to be explained solely 
by the velocity distribution of fragments from a catastrophic collison, given that the parent 
body's diameter could hardly be larger than 50 km. They conclude that, if the Gaulish cluster 
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is indeed coUisional, some other process must have dispersed the fragments after the original 
disruption. It is tempting to suggest that sweeping by the Great-Inequahty resonance could 
have additionally dispersed this group, but we do not think that such a scenario is likely, 
at least in its simplest form. The semimajor axes of the Gaulish cluster members are much 
further apart in velocity terms than their e and i. The Great-Inequality resonance does not 
affect a significantly, but rather tends to induce wide variations in e. In our opinion, the 
present distribution of the members of the Gauhsh cluster could be better explained through 
a depletion by the Great-Inequality resonance, rather than by dispersion. We suggest that 
the Gaulish cluster may have originated from the disruption of a body much larger than 
50 km in diameter. Subsequent sweeping by the Great-Inequality resonance through the 
Gaulish cluster eliminated much of its material. In this view, the three known satellites were 
among those that managed to survive this process, by leaving the resonance before they 
were lost to escape or to collisions with major moons. The erratic nature of any planetary 
migration would offer many possibilities for escape and capture events. The very similar tu 
periods of Erriapo and Tarvos (682 and 674 yrs, respectively) indicate that they would have 
escaped the resonance during the same epoch. 

Direct numerical simulation of the possible effects of planetary migration on the Gaulish 
cluster is clearly needed. In a recent study, Nesvorny et al. (2004) suggest that the high flux 
of impactors in the early Solar System probably dominated the coUisional history of many 
irregular moons. Therefore, it is not completely unexpected that the disruption of the parent 
body might have preceded much of the planetary migration. It is possible that a detailed 
study of the Gauhsh cluster's history could directly constrain the timing of the breakup 
event (if there were one) in relation to the planetary migration. However, while the Great- 
Inequality resonance is very interesting dynamically, its significance is limited to just a few 
of the Saturnian irregulars. At this point, the only conclusion we can make about the history 
of the Gaulish cluster is that its present characteristics arc most likely not primordial, but 
modified by the Great-Inequality resonance. Because of all this, it would be risky to infer 
the origin of the cluster's progenitor based on the present parameters of its members. 



6. Implications for the Origin of Irregular Moons 

Most researchers agree that the irregular satellites are captured bodies, which formed 
in the protosolar nebula independently of the planet. Several capture mechanisms have been 
proposed; the most prominent ones are collisions (Colombo & Franklin 1971), increase in 
planetary mass (Heppenheimer & Porco 1977; Vieira Neto et al. 2004) and aerodynamic drag 
(Pollack et al. 1979). 
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Various aspects of the possible capture and subsequent orbital evolution of Jupiter's 
largest irregular satellite, Himalia, starting from the gas-drag hypothesis, have been explored 
by Cuk & Burns (2004). They find that the capture of Himalia by aerodynamic drag is 
possible, although their model for the Jovian nebula differs from that of Pollack et al. (1979). 
Himalia is arguably the logical starting point if one is interested in exploring the origin of the 
irregular satellites. Since it is both relatively large and prograde, it required more dissipation 
to be permanently captured than any other irregular (with the possible exception of Triton, 
whose origin is outside the scope of this work). If Himalia were captured by dissipation in 
the nebula, such an origin would be even more likely for other irregulars. A more subtle 
difference between Himalia and most other irregulars is that Himalia appears to be outside of 
the region where resonances are common. Cuk & Burns (2004) find that this was likely true 
in the past, since no detectable resonant events can be seen during numerical simulations of 
its post-capture orbital evolution. Using results from previous sections, we can confirm this 
result and put it in the wider context of irregular-satellite dynamics. Fig. 12 is essentially 
the same as Fig. 5, only that now all satellite groups are labeled. The clusters containing 
bodies with diameters larger than 100 km have their names marked with asterisks. It is clear 
that Himalia (in the middle right) lies some distance away from the secular resonance. Not 
only is this true at the present epoch, but Cuk & Burns (2004) argue that the post-capture 
elements of Himalia were probably close to m = 0.13, i = 40°, which is also comfortably 
below the secular resonance. Therefore the result that Himalia's orbital evolution avoided 
resonant passages should not be surprising. 

Globally, Fig. 12 shows significant clustering of the irregular sateUite groups around 
the locations where the secular resonance is possible. We will refer to this super-family as 
the "Main Sequence", due to its superficial visual similarity to the famous feature on the 
Hertzsprung- Russell diagram. This grouping naturally includes all objects that are in (or 
very close to) secular resonance, including Pasiphae (labelled Ps in Fig. 12), Sinope (Sn), 
Siarnaq (Sr), Stephano (St), possibly S/2003 SI and S/2003 J2 (off scale, at m = 0.225), 
as well as some small Jovian satellites that are likely coUisional fragments of Pasiphae and 
Sinope. Almost all Kozai librators are on the "Main Sequence" , too, since the Kozai res- 
onance's boundary (Figs. 4a and 4b) hes usually only several degrees above the secular 
resonance, and the known librators are mostly found just above the boundary. These in- 
clude Kiviuq and Ijiraq (K), Euporie (E), S/2003 J20, S/2003 U3 and possibly S/2002 N2 
and N4. Another dynamical class of objects comprising the "Main Sequence" are reverse- 
circulators. These are the objects whose orbital elements put them between the secular and 
Kozai resonances. Their arguments of pericenter are circulating, but more slowly than their 
nodes. This leads to the precession of w being dominated by 17, which is always in the oppo- 
site direction from the orbital motion (Eq. 5). The satellites exhibiting reverse-circulating 
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behavior are Ananke with its family (A), Themisto (T), Paaliaq (PI), Caliban (Cb), and, 
according to the current orbital solutions, S/2002 Nl and S/2003 Nl. Finally, several ob- 
jects seem to lie relatively close to the secular resonance but still have their w circulating 
regularly in the direction of their orbital motion. At this point, only retrograde objects are 
known to do this, and they include Carme with its cluster, Skadi and probably S/2001 U3. 
Here we need to caution that the orbits of some of the objects we listed above as prospective 
members of the "Main Sequence" , especially the new Neptunian satellites, are still somewhat 
uncertain. However, we believe that it is highly unlikely that there is a systematic error in 
the present orbit solutions that can make them appear closer to the resonances than they 
really are. If anything, a preliminary orbit is likely to miss the resonance. The correction 
to the orbit of S/2003 SI based on observations early in 2004 (Shcppard et al. 2004) puts it 
closer to the secular resonance than the solution plotted in Figs. 5 and 12. 

Table 2 lists the largest members of each suspected irregular sateUite cluster, and reca- 
pitulates our classification of them into dynamical groups. For all objects, the v = —w/Vt 
ratio is also given, in order to demonstrate that our classification has a direct quantitative 
basis. The definition v was chosen so that in the ideal case ofm<< l,e<< 1 and sini << 1, 
V = 1. It can be higher for low-i prograde moons (e.g., for the Moon, v = 2.09), while it is 
generally lower for inclined prograde, and all retrograde, orbits. It is obvious that u = —1 
for Kozai resonance and u — ior the secular resonance. Reverse-circulators have negative 
i/'s (usually close to zero, though) while the objects with zu precessing in the direction of 
orbital motion have u > 1. Among the latter objects, the distinction between those that 
are, and are not, close to the "Main Sequence" is not sharply visible in their u. However, it 
is uniquely defined, with objects having u < 0.2 appearing close to the "Main Sequence" in 
Fig. 12. On the other hand, the classification of irregulars into clusters and the choice as to 
which objects are more-or-less primordial are open to debate (Nesvorny et al. 2003; Grav et 
al. 2003). We decided to be conservative about cluster affiliation: all the irregulars whose 
cluster membership is not clear are taken to be independent objects (i.e., we chose to "split" 
clusters in ambiguous cases). This principle is somewhat softened for the Uranian irregulars, 
where we identify only four "clusters": S/2001 U3 and S/2003 U3 individually, Caliban and 
Stephano together, and all other objects belong to Sycorax's family. Notice that Sycorax's 
escape velocity (which determines the dispersion of fragments) is a much higher fraction of 
its orbital speed than the same ratio for any other irregular, making a very extended cluster 
not too surprising. The extent of the Sycorax cluster was likely increased by post-breakup 
gravitational scattering among family members (Christou 2004). In the case of Neptune, 
by assigning each satellite to its own group, we have most likely overestimated the satellite 
diversity, and it is very probable that the two pairs with similar mean motions might end 
up being genetically related. However, at this time, we believe that such a conclusion would 
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be premature. 

It is interesting to note that none of the three more massive objects (Himaha, Phoebe 
and Sycorax) have orbits anywhere close to the secular resonance. While it is hard to 
claim any statistical significance with only three points, we aver that this dichotomy does 
imply some kind of significant difference that sets the largest few irregulars apart from all 
others. This dichotomy is not surprising because the aerodynamic accelerations - vital in 
satellite capture - act differently on different-sized bodies. In our opinion, large bodies 
were likely captured before the rest, at an epoch when the gas density was high enough to 
cause permanent capture of a 200-km planctcsimal. Cuk & Burns (2004) find no mechanism 
that can arrest the decay of proto-Himalia into Jupiter except for a fast (lO'^-yr timescale) 
evolution of the nebula itself. Such a capture scenario would not result in any preferential 
final orbit for the satellite. 

On the other hand, smaller satellites could not have been present at the same epoch, 
since the strong gas drag would make them rapidly spiral into Jupiter. Therefore, we think 
that the capture of most of the members of the "Main Sequence" postdates that of Himaha, 
Phoebe and Sycorax. The fact that a large fraction of the sateUite groups are found close 
to the secular resonance indicates that some kind of discriminating process must have acted 
on these objects, causing them to end their evolution as members of the "Main Sequence". 
The large amount of orbital evolution that some of these objects experienced (Themisto, 
for example, is very tightly bound) and the existence of several objects that arc exactly in 
the secular resonance, suggests that gas drag played a role in the capture of the smaller 
irregulars, too. Their clustering in the region of slow w precession hints that some process 
acting exclusively on slowly-prcccssing orbits also had a role perhaps in their capture and 
certainly in their evolution. 

Two schemes might explain the "Main Sequence": it might be a region where the satel- 
lites halted their orbital decay into the planet, which is otherwise an unavoidable consequence 
of a sateUite capture in a long-hved gas disk. On the other hand, the "Main Sequence" might 
result from a bias: permanent capture might have been possible but rare at some time in 
the past, with the planetesimals being captured only into certain orbits, which led to the 
clustering of resulting irregular satellites. We are unable unambiguously to endorse either 
of these hypotheses. Nonetheless we discuss each hypothesis briefiy since we have tried 
(inconclusively) to confirm some of their predictions. 

The simplest way of producing the "Main Sequence" would be if the resonant interaction 
supplied energy to objects that otherwise would have decaying orbits. However, such a 
mechanism has to be rejected on theoretical grounds. The argument of secular resonance 
has no dependence on the satellite's mean longitude (see Sec. 4 below Eq. 45; by definition. 
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secular terms cannot include A which varies rapidly). Since the mean longitude (or the 
anomaly, depending on one's choice of Hamiltonian variables) is the conjugate variable to 
the total energy of an orbiting object (Murray & Dermott 1999), only a term containing the 
instantaneous orbital position can be expected to induce changes to an orbit's semimajor 
axis (which is equivalent to the total energy). We are not aware of any published numerical 
experiments that support the contrary conclusion. We also recall here that our preliminary 
numerical experiments on gas-drag evolution (cf. Cuk & Burns 2004) with secular resonances 
never found truly permanent capture. 

Another, subtler way that passage through the secular resonance can stop orbital decay 
is by changing the body's eccentricity. The survival of retrograde satelites, especially Phoebe 
with its low eccentricity orbit, strongly suggests that the circumplanetary disk must have had 
a sharp drop in surface density, most likely near Phoebe's present orbit (for details, see Cuk & 
Burns 2004). In this view, just the pericenter of most sateUites penetrated the nebula so that 
aerodynamic drag affected them only during close approaches. Therefore, a sharp change in 
the eccentricity of the satellite (which would in turn modify the pericenter distance) could 
radically alter the rate of satellite's evolution. If a satellite experienced a significant decrease 
in eccentricity due to secular-resonance passage, it could theoretically become decoupled 
from the disk altogether and end its orbital evolution in the resonance (or very close to it). 
Since Sinope is not only affected by the strongest resonant Hamiltonian (Table 1), but also 
exhibits passages from libration to circulation, we will use this irregular as a test-case for 
this hypothesis. Saha & Tremaine (1993, Fig. 2) have previously described this transition 
between circulation and resonance, and show the results of several integrations from 10^ yr 
to 2 X 10^ yr with all giant planets and with Jupiter alone. Fig. 13a shows the behavior of 
Sinope's eccentricity and its resonant argument (during the simulation shown in Fig. 6b) in 
polar coordinates, i.e., e is shown by the distance from the center, while ^ is represented by 
the angle from the x-axis. Clearly the allowed positions of (e, \t') in this plot lie on two rings 
that overlap at \I' = 0. Motion around the rings corresponds to the two directions of \I''s 
circulation, while the librations are equivalent to crescent-shaped arcs in which the particle 
shifts from one ring to the other without crossing the positive x-axis. Passage through the 
resonance would be seen in this plot as switching from circulation along one ring to the 
opposite motion along the other one. To identify the direction of ^'s circulation along the 
rings. Fig. 13b plots how e changes with time during the same simulation. A comparison 
of Figs. 6b and 13b reveals that the periods of retrograde motion of \1' are correlated with 
high eccentricity, therefore on the outer ring ^ < while \E' > on the inner ring. We 
know that orbits with inclinations below that required for the secular resonance have w 
circulating in the orbital direction, while the higher-inclination ones are reverse-circulators. 
However, since the secular resonance appears to exhibit a downward turn as m increases for 
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retrograde orbits (Fig. 5), the reverse-circulators {vj > 0) are to Sinope's left and objects 
with tz7 < are to the right. Hence, if an object is evolving due to gas drag from larger to 
smaller m (i.e., d < 0), it would increase its eccentricity (and correspondingly accelerate its 
orbital evolution) as it crosses the secular resonance. Therefore we do not believe that this 
mechanism could have arrested the orbital decay of the retrograde satellites. 

For prograde bodies, the shape of the secular resonance (Fig. 5) suggests that a decay- 
ing satellite would cross the resonance from the region of direct precession to that of tu's 
reverse circulation. Fig. 19 in Nesvorny et al. (2003) shows that the instances of retrograde 
circulation of Siarnaq's ^ arc correlated with an eccentricity lower than that seen during 
libration. Therefore, among the prograde satellites, Ae during the resonance passage would 
have been negative, unlike for distant retrograde satellites. However, Fig. 19 of Nesvorny 
et al. (2003) also demonstrates that this change is rather small for Siarnaq, making any 
impact of the resonant passage much less dramatic than it would be for a Sinope-type orbit. 
With this in mind, we conclude that the change in e during the passage through the secular 
resonance does not make a very promising candidate for a process leading to the clustering 
of irregular satellites into the "Main Sequence" . 

The opposite view, that the "Main Sequence" might mirror a bias in the capture mech- 
anism, has some support in the findings of Cuk & Burns (2004). Their Fig. 7 shows that 
the timescalc for temporary capture is related to the precession period of — 'cOpianet- This 
is because capture preferentially happens when the planet is at pericenter (and zero velocity 
curves open up the most; Murray & Dermott 1999, Hamilton & Burns 1992), accordingly 
it results in a temporary satellite orbit having its apocenter close to the Lagrangian point. 
It is easy to see that the combination of these two constraints requires vo — zupianet = at 
the capture. For the duration of temporary capture, the planetesimal experiences resistance 
from the circumplanetary disk during each pericenter passage, and loses orbital energy and 
momentum. If there could be some preferred orbits where temporary capture could last sig- 
nificantly longer than average (usually it persists for about 100 years in the case of Jupiter; 
see the arguments of Cuk & Burns 2004), such objects could become captured even if the 
nebular density were relatively low. In such a nebula, any subsequent collapse of satellite 
orbits due to continuing gas drag would be much slower, and the chance of them surviving 
to the present day would be much higher. 

According to Fig. 12, the region of secular resonances widens towards larger m's for both 
prograde and retrograde orbits, actually becoming rather large for very distant retrograde 
objects. At very large m, satellite orbits become unstable. From the considerations outlined 
immediately above, we would expect that locales in the orbital-element space where the 
"Main Sequence" reaches the stabihty hmit for distant orbits will be places where temporary 
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capture is most likely to become permanent even in the presence of only moderate gas drag. 
In this scenario, a chaotically changing orbit of a temporarily captured object could drift 
into one of these "keyholes" at the ends of the "Main Sequence" . A temporary-satellite orbit 
whose orbital elements put it close to the "Main Sequence" would have a relatively slow 
precession of and this would delay its escape, since w — Wpianet = must be satisfied for 
an orbit to slip out of the Hill sphere's opening. This delay might prolong the temporary 
capture phase by more than an order of magnitude (Table 2), as |z/| < 0.1 is typical for objects 
close to "the Main Sequence" (but not in Kozai resonance) while i/ ~ 1 for other satellites 
(Table 2). After the planetesimal has been permanently captured, its orbital evolution would 
decrease its a and e, while the inclination would lessen for prograde orbits and rise up (i.e., 
move away from i = 180° and toward i = 90°; Cuk & Burns 2004) for retrograde ones. Such 
a migration would force a prograde satellite toward the left and down in Fig. 12, while a 
retrograde one would move to the right and up. This way, the bodies that start at the ends 
of the "Main Sequence" would stay close to it even after substantial orbital migration. In 
such a scheme, all the moons would be captured into orbits with circulating uj but they 
could subsequently evolve into orbits with librating cu. Finally, today's extant objects in 
the secular resonance might have been captured into it later either by very weak residual 
gas drag, sweeping of resonances due to planetary migration (Nesvorny et al. 2003), small 
growth in the planet's mass, or even a variant of the Yarkovsky effect (Cuk et al. 2002b). A 
large population of bodies in the general vicinity of the secular resonance would make such a 
capture likely, even if the changes in the position of the satellites and the resonance relative 
to each other in the orbital-element space was not dramatic. 

Studies of the short-term stability of distant asteroidal satellites by Hamilton & Burns 
(1991) appear to support the "keyhole" hypothesis. They integrate numerous test particles 
starting on circular orbits that exhibit a range of a and i. Their Fig. 15 illustrates the 
initial semimajor axes and inclinations of the particles that survive in orbit for five asteroid 
years (equivalent to about 60 years for the comparable problem with Jupiter). This plot 
clearly shows that there are two inclination values for which stable orbital distances reach 
local maxima: 50°-60° {amax — ^/'^Rh) and 160°-180° {amax — Rh)- We point out that, 
in the context of the Kozai cycle, these are to be considered maximum inclinations, since 
they coincide with minimum eccentricities. Having this in mind, these two distant stabil- 
ity regions are clearly correlated in inclination with our proposed "keyholes". The critical 
distances found by Hamilton & Burns (1991) are much larger than those corresponding to 
long-term stability, but it is likely that the dependence of critical distance on i should be 
somewhat similar. More recently, Vieira Neto & Winter (2001) have studied the stability 
of temporarily captured satellites of Uranus, and their results lend further support to our 
"keyhole" hypothesis. Vieira Neto & Winter (2001) numerically followed orbits of arbitrary 
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eccentricity and inclination, and found a dependence of the capture duration on inclination 
(their Fig. 4b) that is very similar to the results of Hamilton & Burns (1991). Although 
dealing with orbits around objects of very different masses, these numerical experiments 
show that the stability of distant orbits clearly depends on their inclination, and it is likely 
that this dependence will persist for the orbits of permanently captured objects. 

At the present time, we favor the "keyhole" hypothesis as the explanation for the clus- 
tering of small irregular satellites around the secular resonance. Only large-scale direct 
numerical simulations of the capture process could help us decide between the different sce- 
narios. Since we have every reason to think that the likelihood of satellite capture and 
subsequent survival is low, many thousands, if not millions, of particles might need to be 
followed before a statistically significant sample of hypothetical irregular satellites could be 
acquired. The uncertainty in the details of both planetary formation and migration would 
require multiple simulations in order to explore different choices for uncertain parameters. 
Finally, once such hypothetical irregular satellite systems are generated, bodies on inclined 
orbits that are unstable on long timescales (Carruba et al. 2002) would need to be eliminated, 
since their presence would make it harder to compare the synthetic satellite system to the 
natural ones (Astakhov et al. 2003). Such a project requires substantial numerical resources 
and lies outside this paper's scope, which deals primarily with an analytical description of 
today's irregular satellite orbits. 

Apart from the large irregulars, Saturn's "Gaulish" (Albiorix, Erriapo and Tavros) and 
"Norse" (Mundilfari, Suttung, Thrym and Ymir) clusters are also not in the region adjacent 
to the secular resonance. In Sec. 5 we showed that the Gaulish cluster was likely affected 
in the past by the Great-Inequality resonance, which is specific to Saturn, and therefore 
cannot be considered to indicate the primordial irregular-satellite population. We suspect 
that some related resonance might have also affected the Norse cluster, which is not only 
more dispersed in terms of their relative velocities that one would expect for a coUisional 
family, but also exhibits a baffling distribution of eccentricities, which vary more than a 
and i among the cluster members. Clearly, more research is needed to better appreciate how 
planetary migration and the resulting shift in solar-system secular frequencies affect irregular 
satellites, especially those of Saturn (cf. Carruba et al. 2004). One should also recognize 
that other resonances of a similar type, involving either the Great Inequality or perhaps the 
"Lesser Inequality" of Uranus and Neptune (associated with 2Xn — Xu argument), could 
exist at other planets. Alternatively, satellite-satellite scatterings might also be important 
(Christou 2004), especially for members of the "Norse" cluster that share their orbital space 
with Phoebe. 



-36- 



7. Summciry 

In the previous sections, we have shown how a mixed analytical-numerical approach 
can greatly improve our understanding of the irregular-satellite dynamics. Of course, plots 
similar to Figs. 4b and 5 can as well be generated simply by direct numerical integration, and 
doubtless such plots would be more faithful to the behaviors of the real sateUites. However, 
it is noteworthy that the generation of the grid of test particles on which Figs. 4b and 
5 are based, using our secular model, took only about an hour of computing time (on 
an ordinary personal computer with a 2-GHz processor). The same calculation using a 
symplectic integrator would require either several weeks of computing time or a much more 
powerful processor. While such a project would be by no means exceptional, we maintain 
that a reasonably accurate approximation always has certain benefits over a "brute-force" 
approach. Seeing how the Kozai, evection and octupole terms in the disturbing function 
interact and produce the observed distribution of resonances in the orbital-element space 
provides deeper insight into the irregular-satellite dynamics than we would have just by 
finding the resonant locations numerically. 

The dynamics of the irregular satellites has proven to be both exceptionally rich and 
difficult, stimulating a steady flow of publications from several research groups during the 
last few years. This paper deals with many areas of irregular-satellite research; its main 
conclusions are: 

1. The evection and other short-period terms, which do not vanish after averaging over 
the planetary orbital motion, must be included in the disturbing function when constructing 
an accurate secular theory of irregular-satellite motion. 

2. Our model has purely secular and evection terms derived directly, while the total 
contribution caused by other short-period terms is synthesized on the basis of classical lunar 
theories (Tisserand 1894; Brown 1896). This model is found to predict rather accurately 
the precession rate of the longitude of pericenter (uj) for almost all orbital elements. The 
locations of secular resonances can be calculated correctly using this model, while the critical 
inclination for Kozai resonance can be computed for all but very eccentric, retrograde orbits 
(Figs. 4 and 5). 

3. By analyzing the octupole term in the disturbing function (Yokoyama et al. 2003), 
we show that both the locking strength and the location of the stable point of the secular 
resonance can be estimated from relatively simple theoretical considerations, and that these 
predictions agree well with direct numerical integrations of the known irregular satellites' 
motions. 

4. We report a serendipitous discovery of a new secular resonance in the Saturnian 
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system, which consists of a 1:2 commensurabihty between a sateUite's apsidal precession 
frequency with the Great Inequahty of Jupiter and Saturn. This resonance is very strong and 
can make the orbits of satelhtes unstable on short timescales; its location depends sensitively 
on the present orbits of Jupiter and Saturn, and it likely swept through the "Gaulish" cluster 
of Saturn's irregular satellites during planetary migration (cf. Carruba et al. 2004). 

5. We observe that the large majority of irregular-satellite clusters inhabit the region 
close to the secular resonance. We refer to this assemblage as the "Main Sequence" . Only the 
largest satellites (i? > 100 km) and those Saturnians likely affected by the Great-Inequality 
resonance (see preceding point) clearly do not belong to the "Main Sequence" . 

6. We propose a new variant of capture under aerodynamic drag as the mechanism for 
the formation of the "Main Sequence" . Cuk & Burns (2004) have shown that satellite capture 
and escape require an approximate alignment of the apses of planet's and the satellite's orbits. 
Based on this, we argue that temporary capture should last longer for any orbits having slow 
apsidal precession, making them more likely candidates for permanent capture. 

The authors wish to thank Valcrio Carruba, Bob Jacobson, David Nesvorny and Phil 
Nicholson for their help on our irregular satellite dynamics project over the last three years. 
We also thank the anonymous reviewer and especially PDN for comments on an earlier draft. 
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-0.002 I 1 1 1 1 1 1 1 1 1 1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

eccentricity 

Fig. 1. — The dependence of the apsidal precession rate w on sateUite's eccentricity, for a 
near-planar prograde satelhte (m = 0.08, i = 5°) moving about a low-eccentricity pseudo- 
Jupiter (P' = 11.86 yr, but e' = 0.01). Crosses plot the result of a symplectic integration, 
while lines plot the predictions of secular theories. The precession rates for the latter were 
obtained by advancing the secular equations of motion using a Burlisch-Stoer numerical 
integrator. The short-dashed line plots the predictions of Kozai's theory, the long-dashed 
curve also includes the evection terms, while the solid line depicts the results of the mixed 
model derived in this paper (Eqs. 35 and 36). 
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Fig. 2b.— 
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Fig. 2c. — Similar to Fig. 1, for various eccentricities, comparing the motions of both the 
hne of apsides and the hne of nodes predicted by secular theories with direct numerical 
integrations. In all three panels, the precession rate vj for the line of apsides predicted by 
Eq. 35 is plotted as a solid line; likewise, the predicted precession rate for the line of nodes 
(Eq. 36) is plotted by a dotted line. The same two quantities, predicted on the basis of the 
Kozai theory (Eqs. 4 and 5) arc plotted by dashed and dashcd-dottcd lines, respectively. 
All orbits shown have low inclination [i = 5°) and scaled orbital periods of a) m = 0.08 
(prograde), h) m — 0.08 (retrograde), and c) m — 0.16 (retrograde). The planet's period is 
taken to be that of Jupiter (11.86 yrs). 
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-0.01 I 1 1 1 1 1 — ^^^^^ 1 

120 130 140 150 160 170 180 

inclination, degrees 

Fig. 3c. — Similar to Fig. 2, but versus inclination, comparing the motion of both the line of 
apsides and the line of nodes predicted by secular theories with the numerical integrations. 
The continuous lines and the discrete points have the same meanings as in Figs. 2a-c. All 
orbits shown are low eccentricity {emax — 0.2) and have scaled orbital periods of a) m = 0.08 
(prograde), h) m — 0.08 (retrograde), and c) m — 0.16 (retrograde). The minor irregularities 
on some of the curves are artifacts of sampling. 
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Fig. 4a. — The critical inclination (as a function of m) at which Kozai resonance occurs for 
nearly circular orbits, obtained by twenty iterations of Eq. 39 (sohd hne). The individual 
points are the known irregular satellites: Jovians (pluses), Saturnians (crosses), Uranians 
(asterisks) and Neptunians (boxes). The large triangles mark the known Kozai librators, 
according to numerical simulation, whose names are also labeled. 
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m (orbital period, planet=1) 

Fig. 4b. — Critical inclination for the Kozai resonance as determined by numerical advance- 
ment of the secular equations of motion, for particles with a maximum eccentricity of 0.2 
(solid), 0.4 (long dash) and 0.6 (short dash). The symbols for the known individual satellites 
have the same meaning as in panel a). 
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Fig. 5. — The location of the secular resonance determined by numerical advancement of 
the secular equations of motion, for particles with a maximum eccentricity of 0.2 (solid), 0.4 
(long dash) and 0.6 (short dash). The individual points arc the known irregular satellites: 
Jovians (pluses), Saturnians (crosses), Uranians (asterisks) and Ncptunians (boxes). The 
large triangles mark the objects known to be in, or close to, secular resonance, with their 
names also labeled. 
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Fig. 6d. — Evolution of the resonant argument ^ = zu—Wpianet (in radians) during a 300,000- 
yr numerical integration for four resonant or near-resonant satellites: a) Jupiter's Pasiphae, 
b) Jupiter's Sinope, c) Saturn's Siarnaq and d) Uranus's Stephano. 
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Fig. 7b. — Correlation between the resonant argument ^ and the parent planet's eccentricity 
for a) Jupiter's Sinope and b) Saturn's Siarnaq. Note that Sinopc's avoids passing through 
\1/ = when Jupiter's eccentricity is in its upper range, whereas Siarnaq's ^ passes through 
— 77 only when Saturn's eccentricity is at minimum. 
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Fig. 8. — The dependence of S{cosi) as a function of the orbital inchnation (defined here as 
being between 0° and 90°). The values of S'(cosi) and i for four resonant and near-resonant 
objects are also labeled. S{cosi) passes through zero at an average inclination of about 41° 
(139° for retrograde orbits). 
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Fig. 9b. — a) Evolution of eccentricity over a 300,000-yr numerical integration for a particle 
orbiting Saturn whose orbit is affected by the Great-Inequality resonance, b) Evolution of 
the resonant argument ^ = 5A5 — 2\j — 2w — wg during the same integration. 
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Fig. 10. — The dependence of eccentricity variation parameter r] (defined in Eq. 49) on tfie 
precession rate of a particle's uj for 440 test particles integrated numerically over 30,000 yrs 
(see text for details). Secular and the Great Inequality resonances are prominent as regions 
of high rj; the onset of Kozai resonance is also visible, with all the particles to the left of the 
discontinuity being Kozai librators. 
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Fig. 11. — The position of the Great-Inequahty resonance in m — imin space. The top three 
continuous hues plot (in descending order) the Great Inequahty resonance for maximum 
eccentricities of 0.2, 0.4 and 0.6, respectively. The fourth line plots the position of the Great 
Inequality resonance for emax = 0.6, in a planetary configuration that has SA^ — 2Aj = 0.002 
yr~^ (i.e., the Great Inequahty has a period of 500 yrs). The orbital elements of the known 
prograde Saturnian satellites are also shown as individual points. 
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Fig. 12. — Same as Fig. 5, only with all satellite clusters and isolated satellites labeled. 
The groups whose names are followed by an asterisk contain objects with R > 100 km. The 
clustering of groups around the secular resonance is clearly visible. See text for the full 
explanation of the abbrevations used in this figure. 



-56- 



0.4 




-0.4 I 1 1 1 1 1 1 

-0.6 -0.4 -0.2 0.2 0.4 0.6 

e cos(psi) 

Fig. 13a.— 



0.34 




0.2 I 1 1 1 1 1 1 

50000 100000 150000 200000 250000 300000 

time, years 

Fig. 13b. — a) The evolution of Sinope's e and its resonant argument ^ = w — zuj during the 
numerical integration shown in Fig. 6b, plotted in h = ecos\l' and k = esin\l' coordinates. 
The two rings represent two modes of circulation that are characteristic of the regions on 
either side of the secular resonance, while the librations show as crescent-shaped figures that 
include portions of each ring, b) Evolution of e over time in the same simulation. Note that 
the episodes of high-e in this plot are correlated with the instances of retrograde circulation 
of ^ in Fig 6b. 
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Tablc 1: The strength of the resonant term for irregular satelhtes known to be in, or close 
to, secular resonance. 



Satellite 


a/ a' 


< e > 


6i(e) 


< i > 


o 




S'(cosi) 




Pasiphae 


0.03014 


0.4106 


-1.1563 


148.24 


0.11617 


1.0000 


Sinope 


0.03041 


0.2891 


-0.7681 


157.39 


0.22990 


1.3262 


Siarnaq 


0.01251 


0.3041 


-0.8130 


47.81 


-0.07441 


-0.1870 


Stephano 


0.00277 


0.2320 


-0.6033 


143.55 


0.05582 


0.0230 



Table 2: The ratio u = —vj/Vl and our dynamical classification for some irregular satellites. 
Categories are secular resonators (SR), Kozai librators (K), reverse-circulators (RC), "Main 
Sequence" circulators (MSG) and non-"Main Sequence" objects (non-MS). v was computed 
based on our direct numerical integrations described in Sec. 3, except for S/2003 SI, where 
an updated orbital integration was based on Sheppard et al. (2004) 



Satellite 




Class 


SateUite 




Class 


Himalia 


1.10 


non-MS 


Kiviuq 


-1.00 


K 


Pasiphae 


0.00 


SR 


Ymir 


0.28 


non-MS 


Sinope 


0.00 


SR 


Skadi 


0.16 


MSC 


Carme 


0.09 


MSC 


S/2003 SI 


0.00* 


SR (?) 


Ananke 


-0.24 


RC 


Sycorax 


0.35 


non-MS 


Callirrhoe 


-0.26 


RC 


Caliban 


-0.26 


RC 


Themisto 


-0.21 


RC 


S/2001 U3 


0.19 


MSC 


Euporie 


-1.00 


K 


S/2003 U3 


-1.00 


K 


S/2003 J20 


-1.00 


K 


S/2002 Nl 


-0.04 


RC (?) 


Phoebe 


0.59 


non-MS 


S/2002 N2 


-1.00 


K 


Siarnaq 


0.00 


SR 


S/2002 N3 


0.33 


non-MS 


Albiorix 


0.41 


non-MS 


S/2002 N4 


-1.00 


K 


Paaliaq 


-0.08 


RC 


S/2003 Nl 


-0.10 


RC 



